Rearrangement of the vortex lattice due to instabilities of vortex flow 
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With increasing applied current we show that the moving vortex lattice changes its structure 
from a triangular one to a set of parallel vortex rows in a pinning free superconductor. This 
effect originates from the change of the shape of the vortex core due to non-equilibrium effects 
(similar to the mechanism of vortex motion instability in the Larkin-Ovchinnikov theory). The 
moving vortex creates a deficit of quasiparticles in front of its motion and an excess of quasiparticles 
behind the core of the moving vortex. This results in the appearance of a wake (region with 
suppressed order parameter) behind the vortex which attracts other vortices resulting in an effective 
direction-dependent interaction between vortices. When the vortex velocity v reaches the critical 
value v c quasi-phase slip lines (lines with fast vortex motion) appear which may coexist with slowly 
moving vortices between such lines. Our results are found within the framework of the time- 
dependent Ginzburg-Landau equations and are strictly valid when the coherence length £(T) is 
larger or comparable with the decay length Li„ of the non- equilibrium quasiparticle distribution 
function. We qualitatively explain experiments on the instability of vortex flow at low magnetic 
fields when the distance between vortices a ^> Lj n ^> £(T). We speculate that a similar instability 
of the vortex lattice should exist for v > v c even when a < Li„. 

PACS numbers: 74.25.Op, 74.20.De, 73.23.-b 



I. INTRODUCTION 

In 1976 Larkin and Ovchinnikov (LO) [l[ predicted an 
instability of vortex motion that is related to the devi- 
ation of the quasiparticle distribution function from its 
equilibrium value near the vortex core (for review see 
0). When the vortex moves with a velocity v the or- 
der parameter in the vortex core varies on a time scale 
T \ip\ ~ £/ v that can be smaller than the relaxation time 
of the non-equilibrium quasiparticles Ti„ (due to inelastic 
electron-phonon or electron-electron interactions). As a 
consequence, the quasiparticle distribution function de- 
viates from its equilibrium and it results to a shrinkage 
of the vortex core at temperatures close to the critical 
temperature 0]. This effect is mainly connected with 
the removal of quasiparticles from the vortex core by the 
induced electric field 0] and in some respect is similar to 
the dynamic enhancement of superconductivity in weak 
superconducting links 0] • 

Analytical calculations made in the 'dirty' limit (mean 
free path length I of the electrons is smaller than the co- 
herence length £) predicted a decrease of the viscosity co- 
efficient T) of the vortex motion and the existence of a crit- 
ical velocity v c ~ at which the viscous force — vv 
reaches its maximal value 0]- Macroscopically it results 
into a nonlinear current-voltage characteristic V ~ 
with pronounced hysteresis at relatively weak magnetic 
fields [2| • Such a behavior was experimentally observed in 
many low 0,111.13 and hi g h B B EE El El O temper- 
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ature superconductors and quantitative agreement with 
theory was found. From the experimental value of the 
critical voltage V c — v c BL (B is a magnetic flux induc- 
tion, L is the distance between voltage leads) the critical 
velocity v c and relaxation time Ti n were extracted. 

We should stress the nontrivial nature of the LO ef- 
fect. If we use the Bardeen-Stephen expression for 
the viscosity of the vortex motion rj = ^q/2tt^ 2 p n c we 
find that it actually increases if the size of the vortex 
core decreases (with for example a decrease of the tem- 
perature). In the LO theory the vortex core shrinks and 
this results in a decrease of the viscosity coefficient r\. The 
possible explanation of this contradiction is the failure of 
the Bardeen-Stephen model for the vortex as a normal 
cylinder with radius £ in case of a moving vortex with 
high enough velocity. 

At low temperatures small changes in the quasiparti- 
cle distribution function cannot influence the order pa- 
rameter [2|. As a result nonlinear effects will start at 
larger electrical fields and they become significant not 
only in the vortex core but also around the vortex. Effec- 
tively, such non-equilibrium effects were described as due 
to heating of the quasiparticles by the induced electric 
field up to a temperature larger than the sample/phonon 
temperature. This results in a suppression of the order 
parameter near the vortex core and the vortex core ex- 
pands. This effect was used to explain the experimental 
results on the vortex motion instability at low tempera- 
tures for both 'dirty' @, EE El El and 'pure' (I > f) 
superconductors [T{|. In the latter case the instability 
occurs due to a logarithmic dependence of the vortex 
viscosity on the electronic temperature. 

Returning to the LO theory we note that the main as- 
sumption here was the uniform distribution of the non- 
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equilibrium quasiparticle distribution function f(E) in 
the superconductor, (in particular this leads to a field- 
independent critical velocity). From the good quantita- 
tive agreement between theory and experiment in large 
magnetic fields one may conclude that the above con- 
dition is well satisfied when the distance between vor- 
tices satisfies a(B) <C Lin- However experiments at 
low magnetic fields showed that this approach fails when 
a(B) ~ \/®o/B becomes larger than L in @, ©, Qj| 
In Ref. [15| it was proposed that the instability occurs 
when the non-equilibrium distribution becomes uniform 
over the superconductor which is possible if u c Tj„ ~ a(B). 
This leads to a dependence v c ~ yl/B at moderate 
magnetic fields which was observed in many experiments 
0, H, [H, ■ But at lower magnetic fields the critical 
velocity should behave as v c ~ 1/B to explain the field 
independent value for the critical voltage V c [l3[ (the 
same conclusion can be drawn from Fig. 10 of Ref. [3]). 
Note that at these fields the current induced magnetic 
field is still much less than the external magnetic field 
and consequently it cannot explain the observed effect. 

Despite the large number of experimental works there 
is still the fundamental question: what will happen with 
the vortex structure when the critical velocity is ap- 
proached and/or exceeded? In the original paper of 
Larkin and Ovchinnikov it was assumed that the vor- 
tex lattice does not exhibit any structural changes and 
transits for v > v c to a state with a resistance close to the 
normal value (in the current driven regime) . However ex- 
periments on low and high temperature superconductors 
showed that another type of behavior is possible. For 
example, a transition to a state with phase slip lines was 
experimentally observed in Ref. [U SH, [H, HI HI for 
low temperature superconductors and similar IV char- 
acteristics (which were differently interpreted) were ob- 
served in high temperature superconductors in the volt- 
age driven regime [H, [l?], EH- These experimental re- 
sults support the idea that some kind of phase transition 
occurs in the vortex lattice at the instability point and 
regions with fast and slow vortex motion appear in the 
sample 0, M, M M M, 11 SI HI . 

To answer the above questions theoretically one should 
use a rather complicated set of integro-differential equa- 
tions for the order parameter, Green functions of the su- 
perconductor and quasiparticle distribution function 
At the present time it looks almost as an impossible 
task to solve these equations even numerically. There- 
fore, we will limit ourselves to the equations that were 
derived from the microscopic equations for a supercon- 
ductor in the dirty limit under the assumption that the 
relaxation length Li n = \JDTi n of f(E) (D is the diffu- 
sion constant) is smaller than the coherence length £(T) 
[25l [26| . They are the extended (or generalized) time- 
dependent Ginzburg-Landau equations and contain ex- 
plicitly a parameter Tj n . From the very beginning we 
are in a different limit as compared to the LO theory, be- 
cause we consider a non-equilibrium f(E) that is nonuni- 
form in the sample. The longitudinal (odd in energy) 



part f L (E) = f(-E)-f(E) of the non-equilibrium f(E) 
distribution (which is actually responsible for the varia- 
tion of I ■01) is localized only in the region where the time 
derivative d\ip\/dt is finite (see Eqs. (6,10) in Ref. [HI). 
It means that /l(E) is non zero only near the moving 
vortex core. 




FIG. 1: Model system - a superconducting slab (infinitely 
long in the y and z directions) in a parallel magnetic field H 
with transport current /. 



The paper is organized as follows. In section II we 
present our model system. In section III we study 
the rearrangement of the vortex structure due to non- 
equilibrium effects at moderate and high magnetic fields 
and in Sec. IV the non-equilibrium vortex dynamics at 
zero magnetic field. Finally, in section V we discuss our 
results and make a comparison with experiments. 



II. MODEL SYSTEM 

As a model system we use a bulk superconductor which 
is infinite in the z and y directions and is finite in the x- 
direction (see Fig. 1). This model allows us to study 
the effect of the nonuniform current distribution in the 
superconductor (due to transport current) on the vortex 
dynamics at zero and low magnetic fields. 

In our calculations wc neglect the possibility of the for- 
mation of curved vortices in the z direction and therefore 
our problem becomes two-dimensional. The generalized 
time dependent Ginzburg-Landau equations in our case 
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can be written as 
u 



0.05 



d_ 7 2 diif 
V /1 + 7 2 |V'| 2 + ~2~dt~ 

(V - iA)V + (1 - M'M 



dA 
~dt 



= Rc h£>*(-iV - A)V>] - ^ 2 rot rotA, 



(la) 



(lb) 



where the parameter 7 = 2teA(T)/7i is the product 
of the inelastic collision time te for clcctron-phonon 
scattering and A(T) = ^TcW 1 / 2 /^! - T/T c is the 
value of the order parameter at temperature T which 
follows from Gor'kov's derivation [27j of the Ginzburg- 
Landau equations. In Eqs. l(a,b) the physical quanti- 
ties are measured in dimensionless units: temperature 
in units of the critical temperature T c , the vector po- 
tential A = (A x ,A y ,0) and the momentum of the su- 
perconducting condensate p = V</> — A are scaled in 
units <5>o/(27r£(T)) (where 3>o is the quantum of mag- 
netic flux), the order parameter ip — \ip\e 1 ^ in units of 
A(T) and the coordinates are in units of the coherence 
length f(T) = (8k B (T c - Tj/nhD)- 1 / 2 . Time is scaled 
in units of the Ginzburg-Landau relaxation time tgl — 
irh/8kB{T c — T)u, voltage (V) is in units of ipo = h/2ercL 
(a n is the normal-state conductivity). In these units the 
magnetic field is scaled with H c2 = <I>o/27r£ 2 and the cur- 
rent density with jo — <r n h/2eTGLS,- Following Ref. [lf| 
the parameter u is taken to be equal to 5.79. 

Instead of the usual gauge divA = 0we chose the elec- 
trostatic potential equal to zero ip = 0. The semi- implicit 
algorithm was used [28| which provides an effective nu- 
merical solution of Eqs. l(a,b) for the case of large k 
values. We apply periodic boundary conditions in the y 
direction ip(y) = ip(y + L), A(y)=A(y+L) (L is the pe- 
riod - see Fig. 1) and the superconductor- vacuum bound- 
ary conditions in the x direction (V x — iA x )il>\ x =o.w — 0. 
The transport current was introduced via the boundary 
condition for the vector potential in the x - direction 
rotA| z (x = 0, W) = H ± Hi where H I = 2wl/c is the 
magnetic field induced by the current / (per unit length 
in the z-direction) and H is the applied magnetic field. 
In all our calculations we chose k = 5 and the parameter 
7 is varied from to 40. 

Due to the discrete nature of the vortex motion the 
voltage response in our system is a time-dependent vari- 
able. We average it over a finite time interval which is 
taken to be larger than the period of the voltage vari- 
ation. But this time interval can be comparable to the 
switching time between different dynamic phases and it 
smoothens the current-voltage characteristics at the tran- 
sition points. 



III. REARRANGEMENT OF THE VORTEX 
LATTICE 

First we consider the situation when the applied mag- 
netic field is larger than the magnetic field due to the 
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FIG. 2: (Color online) Current-voltage characteristics of the 
superconducting slab with width W = 50£, 7 = 10 and H = 
0.3-ff C 2. Current increases from zero to some finite value. In 
the inset snap-shots of the order parameter at different values 
of the applied current are showed. 



transport current H 3> Hi. Therefore the current den- 
sity distribution is almost uniform over the width of 
the superconductor and it simplifies the analysis of the 
obtained numerical results. In Fig. 2 we present the 
current-voltage (IV) characteristic of our system (W = 
50£, L = 25£, 7 = 10) under investigation at H = Q.3H C 2- 
Vortex flow starts at some finite current / s (due to the 
presence of the surface barrier in the system) and the 
vortex structure is close to the triangular lattice. With 
increasing current it transforms to a row-like structure 
but keeping the triangular ordering (see point 2 in Fig. 
2). Increasing the current (arrow 2 in Fig. 2) there is 
a transition in the vortex structure which is visible as 
a kink in the IV characteristic. The number of vortex 
rows decreases (from 6 to 5 in this particular case) and 
the number of vortices in the rows increases. Note that 
the number of vortices in the system does not change 
and the kink in the IV characteristic is connected with a 
faster vortex motion in this vortex configuration. At the 
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current indicated by the arrow 4 there is a second transi- 
tion where the number of rows decreases further (from 5 
to 4) and the distance between the vortices in each row 
decreases. Simultaneously the vortex velocity increases 
steeply and we have a transition to a state with a much 
larger voltage. 

The transitions in the vortex lattice will be explained 
by the modification of the shape of the vortex core due to 
non-equilibrium processes. Indeed, the motion of the vor- 
tex means a suppression of the order parameter in front 
of the vortex and recovering the order parameter behind 
it (see Fig. 3). If the vortex velocity is large enough 
(v < C/ T in) the number of quasiparticles in front of the 
vortex will be less than the equilibrium value and larger 
behind the vortex due to the finite relaxation time Ti n 
of the quasiparticle distribution function. Effectively, we 
have a cooling of the quasiparticles in front of the vortex 
and heating behind the vortex (see Fig. 3). This effect is 
very similar to the behavior of a superconducting weak 
link at voltages V < 1 /Ti n [29|, [3(| [Hl| when there is a 
cooling at the decrease of the order parameter and a heat- 
ing when the order parameter increases in the weak link. 
Because the relaxation time of the order parameter de- 
pends on the temperature as t\m ~ 1/(T C — T) we have a 
long healing time of the order parameter behind the vor- 
tex and a short time suppression of the order parameter 
in front of the vortex. It leads to an elongated shape of 
the vortex core with a point where = shifted to the 
direction of the vortex motion. This is visible (see Fig. 
2) for vortices close to the right side of the superconduc- 
tor, where the current density and the vortex velocity are 
maximal. 

When the transition from 5 to 4 vortex rows in the 
vortex lattice occurs the distance between the vortices 
suddenly decreases. If the vortex velocity is large enough 
such that v > v c ~ o/tui (a is the distance between 
vortices in the row) the order parameter does not have 
sufficient time to recover after every vortex passage in 
the row and will be strongly suppressed along the 
vortex trajectory. It speeds up the vortex motion because 
the time variation of depends on the value of 
r |^| ~ l\ip\ T GL [32| . This is the reason for a transition 
from slow to fast vortex motion (quasi-phase slip line 
behavior) and a steep increase in the voltage at the point 
where the current is I c and the voltage is V c in Fig. 2. 

At higher magnetic fields the situation is similar to 
the case H = 0.3H C 2 but in addition a transition to a 
state with vortex rows moving with different velocities is 
possible (see Fig. 4). When a fast vortex row (which we 
will further call a quasi-phase slip line (PSL)) appears in 
the sample the superconducting current decreases around 
the PSL on the scale of the decay of the electric field E 
(or charge imbalance) Le Q. Then vortices in adjacent 
to PSL areas are forced to move with smaller velocities 
because the superconducting current mainly drives them. 
In the framework of the model equations l(a,b) Le — 
\/j/ u £, > L in [32| and we found indeed that for larger 
values of 7 the current and magnetic field range over 
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FIG. 3: Deformation of the vortex core due to vortex mo- 
tion (schematic). The density of the horizontal lines shows 
the density of the quasiparticles. In case the diffusion length 
Li n is smaller than the coherence length £(T) diffusion of the 
quasiparticles is not strong and locally there is an effective 
cooling and heating of the quasiparticles. 



which this structure may exist increases. For example, 
for 7 = 20 the slow and fast vortex rows may coexist 
already at H = Q.2>H C 2 and for 7 = 40 they may coexist 
even at H = 0. 

Note, that in contrast to the case H — 0.3H C 2 the in- 
stability of the vortex lattice leads to quasi-chaotic vortex 
motion at currents between points 2 and 3. We relate this 
to the usage of periodic boundary conditions. For exam- 
ple in case of H = Q.3H C 2 the same chaotic behavior (not 
shown here) disappears between points 2 and 3 in Fig. 2 
with an increase of the period of our system by a factor 
of two (with a small effect on the values of the currents 
where structural transitions occur). But for H — QAH C 2 
doubling the period did not result into any effect. 

We explain the influence of the boundary conditions by 
incommensurability effects between the period L and Le- 
Actually the latter length defines the scale of the interac- 
tion between phase slip lines. Changing the parameter 7 




FIG. 4: (Color online) Current-voltage characteristics of the 
superconducting slab with width W = 50£, 7 = 10 and H = 
0AH c2 . 



FIG. 5: (Color online) Current-voltage characteristics of the 
superconducting slab for the parameters of Fig. 2 and H = 
0.7H c2 . 



we change Le- For example for 7 = 20 and H = 0AH C 2 
we did not observe any irregular vortex distribution for 
the superconductor for the parameters corresponding to 
Fig. 3 even when L = 25£. 

At magnetic fields close to H C 2 there are also transi- 
tions in the vortex structures, but they are masked by a 
large number of possible transitions due to the increased 
number of vortices in the system (see Fig. 5). The kinks 
in the current-voltage characteristics become almost in- 
visible and the jumps in the voltage gradually decreases 
at the current I c where the quasi-phase slip lines appear 
in the system. 



IV. VORTEX MOTION AT ZERO MAGNETIC 
FIELD 

In Fig. 6 we present the IV characteristic of the same 
sample as in Fig. 2 at zero magnetic field. At low cur- 
rents we have slow vortex motion while at large current 
quasi-phase slip lines appears. We should note that we 
did not observe any structural changes in the vortex lat- 
tice at low magnetic fields due to the small number of 



vortices and hence the large distance between them. 

In case of zero magnetic field the current density is 
strongly nonuniform over the width of the sample (see 
Fig. 7). When the current exceeds the critical value 
I s (current of suppression of the surface barrier for vor- 
tex entry) the Meissner state is destroyed, vortices and 
antivortices enter the sample, pass through it and annihi- 
late in the center. This process results in the appearance 
of an additional maximum in the current density in the 
center of the sample - see Fig. 7 (in agreement with ana- 
lytical calculations of Refs. [33, [3J]). For sample param- 
eters of Fig. 6 the quasi-phase slip behavior starts when 
the current density in the center reaches the value close 
to the depairing current density jdep- We found that at 
this moment the annihilation of vortex-antivortex pairs 
speeds up and it provides a favourable condition for fast 
vortex motion across the whole sample. However with 
increasing width of the sample the transition to the fast 
vortex motion behavior starts at a larger current (see 
Fig. 8) . As the speed of the fleet is defined by the speed 
of the slowest ship the nucleation of the quasi-phase slip 
line depends on the vortex motion in the place where the 
current density is minimal. When we increase the width 
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FIG. 6: (Color online) Current-voltage characteristics of the 
superconducting slab with the same parameters as in Fig. 2, 
but now for zero magnetic field. 

of the sample we decrease the minimal current density 
jmin in the sample (compare Figs. 8(b) and 7). When 
jmin reaches the critical value the quasi phase slip line 
nucleates in the sample. This critical value is smaller the 
larger 7. This result is closely connected with the find- 
ings of Ref. where it was shown that the phase slip 
process does not exist in a quasi-lD superconductor and 
2D thin superconducting films with uniform distribution 
of the current density [35j if the current density is smaller 
than some critical value j cl (7) . 

We like to stress that we did not find that the vor- 
tices and anti-vortices can pass through each other as 
predicted in Ref. [36| . Probably the uniform current dis- 
tribution used in Ref. [3(| is essential to observe this 
effect. 



V. DISCUSSION 

A. Comparison with other theoretical works 

Our results strongly supp ort the intuitive idea (pub- 
lished already in Ref. [20() about nucleation of phase 



FIG. 7: (Color online) Distribution of the time and length 
averaged current density < j >= j n + j s over the width of 
the superconducting slab at zero magnetic field and different 
values of the transport current. Numbers 1-3 in the figure 
corresponds to different values of the transport current in Fig. 
6. 



slip lines at large currents against a background of vortex 
flow. This idea was further developed in the theoretical 
work [37| where the current-voltage characteristics of a 
wide superconducting film with both viscous vortex flow 
and phase slip lines was calculated. However the author 
used equations that are averaged over the inter-vortex 
distance and did not find the rearrangements of the vor- 
tex structure at high vortex velocities. 

In Ref. [38| the appearance of the wake behind the 
moving vortex was theoretically predicted on the basis of 
an analytical solution of the linearized equation (la) for 
the absolute value of the order parameter. Actually such 
a wake should exist even in the simple time-dependent 
Ginzburg-Landau equation (with 7 = 0) due to the finite 
time for the order parameter relaxation 71^1 ~ tql- In- 
deed, when a vortex moves the current density in front 
of its motion is the sum of the current density from the 
vortex j vor t and the transport current density j ext and 
likewise behind the vortex it is the difference j vor t — jext- 
The time relaxation of the order parameter depends on 
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FIG. 8: (Color online) (a) Current-voltage characteristic of 
the superconducting slab with parameters: 7 = 20, W = 75£ 
at zero magnetic field. In the inset to Fig. (a) we present snap- 
shots of \il>\ at different values of the current, (b) Distribution 
of the time and length averaged current density < j >= j n +js 
for wider sample. Numbers correspond to different values of 
the current at Fig. (a). 

the value of the current density if j is close to the depair- 
ing current density (the larger the current the smaller is 
Tim - see chapter 11.4 in Ref. 4]). If j ext is close to 
jdep the difference between ti ,/, i in front and behind the 
moving vortex is substantial [391 ] and the moving vortex 
becomes elongated in the direction of its motion. 

The change in the shape of the vortex was also found 
in Ref. [4(| on the basis of a numerical solution of the 
2D time-dependent Ginzburg-Landau equations. Such 
vortices were called kinematic vortices due to their high 
velocity. They were found to exist when a quasi-phasc 
slip line nucleated in the sample. The system resembles 
a Josephson vortex in a long Josephson junction where 
anisotropy is connected with different penetration depths 
of the screening current along and across the Joseph- 



son junction (see chapter 6.4 in [4}). In case of a phase 
slip line the anisotropy is connected with a strongly sup- 
pressed order parameter in the direction of the vortex 
motion (along the quasi-phase slip line). 

In both of the above works the dependence of the re- 
laxation time t\^\ on the value of the order parameter 
was ignored. In Ref. [4(| the term 7 2 5|?/>| 2 /2<9i on the 
left hand side of Eq. 1(a) was neglected and instead of 
the coefficient uj 'a/1 + 7 2 |V'| 2 the variable parameter u* 
was used Actually in Ref. [38] the same approach 
was used because in Eq. 1(a) a fixed value for the order 
parameter |V>| in the term T|^,| = tqlu^JI + 7 2 |"0| 2 was 
put. 

But in a TDGL model with 7 = and arbitrary value 
of u we did not find any steep transition from slow vortex 
flow to fast vortex motion (quasi-phase slip line) at finite 
value of the applied magnetic field. The reason is simple: 
in that model the relaxation time of the order parameter 
practically does not depend on the value of the order 
parameter and the mechanism of the switching in vortex 
motion discussed in Sec. Ill does not work. Besides we 
did not find any vortex structure rearrangement in the 
model with 7 = and arbitrary u. Probably the change 
of the shape of the vortex is small in the above simplified 
model. 

In Ref. [35| Eq. 1(a) was coupled with the equation 
for the electrostatic potential and the transition from the 
slow vortex flow to the quasi-phase slip lines behavior was 
numerically observed in case of thin 2D superconducting 
films of finite length in a perpendicular magnetic field. 
However the rearrangement of the vortex lattice and the 
coexistence of the fast and slow vortex motion were not 
found because of the small width of the samples. 

B. Range of validity of the obtained results 

Our results are strictly valid only when Li n < £(T) 
while usually in experiments Li n £(T). But it is obvi- 
ous that cooling and heating of the quasiparticles around 
the vortex core occurs in both limits. Large Li n provides 
some kind of space averaging of these different processes 
(due to diffusion of the non-equilibrium quasiparticles 
from the overheated region to the overcooled one). In 
the framework of the LO theory the zero order effect was 
calculated (when the dependence on the direction of the 
vortex motion was neglected) which roughly leads to an 
effective cooling of the system and a symmetrical shrink- 
age of the vortex core if the distance between the vortices 
is much smaller than Li n . In this limit the separation of 
the system into slow moving vortices and quasi-phase slip 
lines is in principle possible for v > v c when the vortex 
motion becomes unstable. The origin for such a behavior 
is the presence of the normal component of the current 
density (electric field) and its finite decay length from 
the phase slip line Le- The slow vortex motion between 
the quasi-phase slip line may occur due to a weakening of 
the superconducting component of the current near the 
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quasi-phase slip line. In the LO theory this effect was 
neglected and only deviations of the longitudinal part of 
the quasiparticle distribution function from equilibrium 
was taken into account while the transverse (even in en- 
ergy) part of f(E) is responsible for the appearance of 
the finite normal current in the superconductor 

The situation is different if the distance between the 
vortices exceeds Li n and the effective averaging becomes 
weaker. The anisotropy of the vortex core and the ef- 
fective attraction between vortices should be more pro- 
nounced and leads to the appearance of vortex rows/lines 
(slow or fast) even at vortex velocities less than v c . 
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C. Comparison with experiments 



The important property which follows from our cal- 
culations is the weak dependence of the critical voltage 
V = V c on the applied magnetic field (see inset in Fig. 
9). We explained it by the rearrangements of the vortex 
lattice when the vortex velocity approaches v c . Because 
in this case we do not have a triangular vortex lattice 
the distance betwe en vortices in the rows will be smaller 
than a ~ ^JWqJB and defined by a ~ 1/B dependence. 
Indeed, at current I = I c the transition to a state with 
4 quasi-phase slip lines occurs (see Figs. 2,4,5) in the 
magnetic field range 0.3-0.7 H C 2 while the number of the 
vortices in the sample increases linearly with magnetic 
field. Assuming that the transition to the fast vortex 
motion state occurs when the distance between the vor- 
tices in the row is equal to a ~ 1/B ~ v c t\^\ (see Sec. 
Ill) we obtain v c ~ 1/Bt^ and the field independent 
critical voltage V c — v c BL. Such a dependence was ex- 
perimentally observed in Ref. [l3| for both low and high 
temperature superconductors in the low magnetic field 
regime where the vortex separation a 3> Li n and the 
self-field of the transport current was negligible. 

From the above estimations it follows that V c ~ 1/7. 
We see from the inset in Fig. 9 that indeed V c ~ 1/7 
and besides the resistivity of the superconductor at low 
currents follows the dependence p/p n ~ I/7 (see Fig. 
9) analytically found in 41] for large magnetic fields in 
the temperature interval T c — h/ksTin < T < T c where 
equations l(a,b) are valid. 

Our estimation of the critical velocity was made in 
the spirit of the paper of Doettinger et al. [HI]. They 
supposed that the vortex motion instability occurs when 
the non-equilibrium quasiparticles induced in the vortex 
core does not have time to relax to equilibrium when the 
next vortex arrives to the place where they were induced 
Tin = a/v c . Actually it means that the order parameter 
did not have time to increase (because its value depends 
on f{E) and it cannot grow faster than Tj n ). 

The stair-like structure of the IV characteristics (which 
is a fingerprint of the nucleation of the phase slip centers 
or quasi-phase slip lines) was observed both in low and 
high temperature superconductors [2(| HH, [H, [H, [24[ 
at low magnetic fields. In wide samples (in which the 



FIG. 9: (Color online) Current-voltage characteristics of the 
superconducting slab with width W = 50£, k — 5 and differ- 
ent values of the parameter 7 at two values of the magnetic 
field. In the inset the dependence of the the critical voltage on 
the magnetic field is shown for three values of the parameter 
7 = 10,20,40. 



strongly nonuniform current density distribution over the 
width of the sample is realized in the Meissner state) 
a slow vortex motion was found low currents which is 
changed into the quasi-phases slip line behavior at higher 
currents [2(J, Hi], H3, HI1 • It is interesting to note that the 
stair-like structure was also experimentally observed at 
high magnetic fields (and temperatures far below critical 
one) in high-temperature superconductors [13, 13 ■ I n 
this case one was able to observe it only in the voltage 
driven regime. 

In the experiments of Kunchur et. al [l6|, [l7|, EH the 
quasi-phase slip lines become visible because the penetra- 
tion length of the electric field L e increases in high mag- 
netic field due to the suppression of the order parameter 
by vortices and a steep increase of Ti n at low tempera- 
tures @. At low magnetic fields and high temperatures 
too many quasi-phase slip lines appear simultaneously at 
I = I c and it smooths out the stair structure of the IV 
characteristic. 

We explain the absence of the stair-like structure in 
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the majority of the experiments where the LO effect was 
studied due to the large length of the samples as com- 
pared to Le- At the instability point many quasi- phase 
slip lines should appear in such a sample and it results 
into a large heating of the electronic subsystem and the 
sample itself. That could be the reason for the transition 
to the nearly normal state. Besides in those samples the 
role of the voltage can be very pronounced and it leads 
to an additional suppression of superconductivity, the ef- 
fect which is absent in the model equations l(a,b) 
Therefore it would be interesting to perform an experi- 
ment on a short superconducting bridge with length of 
about several Le at different temperatures and magnetic 
fields. Taking into account the strong dependence of r in 
and hence Le on temperature it would be wort to study 
several samples with different length satisfying the con- 
dition L ~ Le at different temperatures. By variation of 
the magnetic field one could observe the predicted split- 
ting of the vortex flow into regions with fast vortex flow 
(quasi-phase slip lines) and slow vortex flow. A good can- 
didate is NbGe which has rather low bulk pinning even 
atT~T c /2 0|. 



D. Hysteretic behavior 

The hysteresis is almost absent in zero magnetic field 
for a superconducting slab with W = 50, k = 5 and 7 
1 (see Fig. 10). Hysteresis appears when we increase 
the applied magnetic field (see Fig. 10) or decrease the 
width of the sample. In both cases the current density 
distribution becomes more uniform in the sample and it 
brings hysteresis in the system. 

The physical reason for this effect is as follows. Con- 
sider at first the case of zero applied magnetic field 
H = 0. In Ref. [H, [H[ it was found that in a super- 
conductor with uniform current density distribution the 
phase slip center/line does not exist at current density 
jdil) which can be smaller than jde-p (in case of zero 
fluctuations). But the superconducting state in such a 
system can be stable up to I c = jdepdW. When a quasi- 
phase slip line nucleates at I c it can exist up to a smaller 
current I\ = j c idW and it leads to hysteretic behavior. 
When we take into account screening effects, the current 
density distribution becomes nonuniform over the width 
of the sample. It is maximal on the edge and minimal 
in the center of the sample being in the Meissner state. 
At current I — I s < I c the current density on the edge 
reaches the depairing current density and the supercon- 
ducting state becomes unstable. Vortices enter the sam- 
ple and if the minimal current density is larger than j c \ 
they move fast and form a quasi-phase slip line. If the 
minimal current density is smaller than j c \ they move 
slowly and form quasi-phase slip lines at a larger current 
when the condition jVnin > jci is fulfilled. In the latter 
case hysteresis is absent because the transition from the 
slow to the fast vortex motion (or vice versus) occurs 
when the current density in one point (over the width of 



the sample) reaches the critical value. It does not lead 
to a crucial redistribution of the normal and supercon- 
ducting current density over the sample and the vortex 
motion is almost non-hysteretic. 

When a magnetic field is applied the transition to 
the fast vortex motion occurs when the current density 
reaches a critical value (which depends on the magnetic 
field) practically over the whole sample. As a result the 
distribution of the normal and superconducting current 
density changes drastically over the whole sample and it 
provides the hysteretic behavior for vortex motion in our 
model system. 




0.00 0.05 0.10 0.15 0.20 0.25 



i/G w) 



FIG. 10: (Color online) Hysteresis in the current- voltage char- 
acteristics of the superconducting slab with W = 50£, 7 = 10 
and different magnetic fields. 



In our calculations we did not take into account heat- 
ing effects. Their incorporation in the considered model 
will increase the hysteresis [43| and mask all effects [44| 
discussed in our paper. 
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E. Synchronization of vortex motion in adjacent 
vortex rows 

Can the motion of the vortices in two adjacent quasi- 
phase slip rows/lines be synchronized? It was found in 
many papers (for a review see (4o1. |46|) that the dynam- 
ics of the order parameter in one phase slip center may 
influence the dynamics of the order parameter in a re- 
mote phase slip center even if the distance between them 
is large. The effect is mainly connected with the long 
decay length of the quasiparticle (normal) current from 
the phase slip center. As a result the ac component of 
the normal current affects the oscillations of the order 
parameter in the other phase slip center in a way similar 
to a Josephson junction under the action of an external 
ac current or electromagnetic radiation. The interaction 
between phase slip centers becomes more complicated if 
one take into account the nucleation of the charge imbal- 
ance waves [45| which can both enhance and suppress the 
synchronization of the the order parameter oscillations in 
adjacent phase slip centers/lines [45|. 




-o.oi 
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FIG. 11: (Color online) Time dependence of the instant volt- 
age and local electric field at different values of the applied 
current for the sample with parameters as in Fig. 2. 

In our calculations, effect due to charge imbalance 
waves were not considered and quasi-phase slip lines in- 
teract only via the ac component of the normal current. 



In Fig. 11 we present the time dependence of the aver- 
aged, over the length of the sample, electric field at the 
edge of the superconductor and the local electric field in 
the point where the quasi-phase slip line at large current 
nucleates (at H = 0.3H C 2 and different currents). At low 
currents, when the induced voltage and normal current 
are small the motion of the vortices in adjacent rows is 
out of phase because of weak interactions between vor- 
tex and strong repulsion between vortices. Contrary, at 
large currents the exit of the vortices from the supercon- 
ductor in adjacent rows becomes in phase. It means that 
the emitted electromagnetic radiation should be consid- 
erably enhanced in this case in comparison with the low 
current limit. 

The frequency of the radiation can be tuned by an 
applied magnetic field or /and by applied current. The 
applied magnetic field changes the number of the vor- 
tices in the row and hence changes the distance between 
vortices while transport current changes the vortex veloc- 
ity. Both factors influences the frequency of the extracted 
radiation v = v/a. Taking typical values for the critical 
velocity in YBCO at T=77 K and B=l T v c = 10 3 m/s 
[§] and the inter- vortex distance a ~ y/& /B — 40nm we 
obtain v ~ 4 • 10 11 Hz. In the quasi-phase slip line regime 
the vortex velocity v^> v c and the frequency approaches 
the THz regime. 



VI. CONCLUSIONS 

In the framework of the generalized time-dependent 
Ginzburg-Landau equations we showed that with increas- 
ing applied current the moving Abrikosov vortex lattice 
changes its structure from a triangular one to a set of par- 
allel vortex rows. The effect originates from changes in 
the shape of the moving vortex. The vortex core becomes 
elongated in the direction of vortex motion because of 
different relaxation times of the order parameter in front 
and behind the moving vortex. In front of the moving 
vortex the order parameter may vary very fast due to a 
large value of the local current density and a deficit of 
quasiparticles in comparison with its equilibrium value. 
Contrary the number of the quasiparticles exceeds locally 
their equilibrium value and the current density is small 
behind the moving vortex and it increases the relaxation 
time of the order parameter. This results in the appear- 
ance of a wake behind the vortex which attracts other 
vortices. 

We found that the rearrangement of the vortex lattice 
results in field-independent value of the critical voltage 
at which the transition to the state with quasi-phase slip 
lines occur. This is connected with changes of the inter- 
vortex distance at the structural transitions of the vortex 
lattice. In a triangular lattice the average distance be- 
tween vortices varies as a ~ 1/ \AB while in case of vortex 
rows the minimal inter- vortex distance decreases with in- 
creasing magnetic field as a ~ 1/B. It results in depen- 
dence V c ~ \/B for triangular lattice and V c ~ const for 
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vortex rows. 

At some magnetic field the quasi-phase slip lines can 
coexist with slowly moving vortices between such lines. 
Besides we found that the motion of the vortices in ad- 
jacent quasi-phase slip lines can be synchronized at large 
vortex velocity v > v c . Both effects are possible due 
to the long decay length of the normal current near the 
quasi-phase slip line. It decreases the superconducting 
component of the current in the system and provides syn- 
chronization of oscillations in the order parameter at the 
quasi-phase slip lines. 

Although our results are strictly valid when £(T) > Li n 
they qualitatively explain experiments on the instability 
of the vortex flow at low magnetic fields when the dis- 



tance between vortices a 3> L in ^> £(T). Besides our 
results support the idea that a similar instability of the 
vortex lattice should exist for v > v c even when a < L- m . 
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